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ABSTRACT 


Repeated  surface  wave  group  velocity  measurements  over  the  same  paths 
were  made  by  using  large  aftershocks  of  several  great  (M  > 7)  earthquakes 
recently  occurred  in  China.  Multiple  filter  technique  was  applied  to  the 
properly  rotated  three-component  digital  data  from  Seismological  Research 
Observatory  (SRO)  stations,  so  that  both  Rayleigh-  and  Love-wave  dispersion 
data  were  obtained  over  a number  of  paths  crossing  various  tectonic  provinces 
of  China.  In  several  cases,  higher  mode  data  were  also  derived.  An  estimate 
of  uncertainty  of  these  dispersion  data  was  obtained  from  the  repeated  measure- 
ments with  identical  source-receiver  geometry.  The  generalized  surface  wave 
inversion  technique  was  applied  to  these  multi-mode  dispersion  data,  and 
crustal  and  upper  mantle  structures  were  derived  for  various  tectonic  provinces 
of  China.  The  results  clearly  demonstrate  that  the  Chinese  Mainland  is  far 
from  being  laterally  homogeneous,  the  lateral  heterogeneities  are  closely 
related  to  the  past  tectonic  developments.  A particularly  unusual  crustal 
and  upper  mantle  structure  is  found  underlying  the  Tibet  plateau. 


(1)  Chinghai -Tibet  (southwestern  China), 

(2)  North  China, 

(3)  South  China,  and 


(4)  Northwestern  China  (including  the  Tien  Shan  fold  belt). 

These  subdivisions  are  similar  to  those  of  McElhinny  (1973),  based  on  paleomagnetic 
data,  Whittington  and  Hughes  (1972)  and  Jell  (1973),  based  on  paleontologic  data, 
and  Dewey  and  Bird  (1970),  Hamilton  (1970)  and  Burrett  (1974),  based  on  regional 
geology.  Each  subplate  has  different  geomorphologic  features,  gravity  anomalies 
and  crustal  thickness  (Tung,  1974;  Sun  and  Teng,  1977).  Major  seismic  events 
occur  along  the  subplate  boundaries.  This  seismicity  is  related  to  the  large- 
scale  tectonic  movements  discussed  by  Molnar  and  Tapponier  (1975,  1 97JZ1 . 

The  division  of  China  into  four  major  subplates  is  helpful  from  the  standpoint 
of  surface  wave  dispersion  studies.  The  concept  of  regionalization,  first  pro- 
posed by  Toksoz  and  Anderson  (1966),  consists  of  dividing  a surface  wave  path 
into  oceanic,  shield  and  mountain  tectonic  regions.  Each  region  has  a different 
effect  upon  the  average  group  velocity.  Toksoz  and  Anderson  (1966)  noted  that 
shield  areas  raise  the  average  group  velocity,  whereas  active  tectonic  regions 
have  the  effect  of  lowering  the  average  group  velocity.  As  a result,  surface  wave 
paths  crossing  different  tectonic  regions  should  show  different  trends  in  their 
average  group  velocity  dispersion  curves.  The  relationship  between  the  concept 
of  regionalization  and  its  tectonic  implications  has  been  studied  (Tung,  1974); 
the  results  show  that  it  is  valid  to  divide  China  tectonically  into  several  sub- 
plates, so  that  pure  path  surface  wave  dispersion  data  can  be  derived  for  crustal 
and  upper  mantle  structure  inversion.  This  paper  makes  use  of  the  digital  data 
that  has  become  available  due  to  the  recent  occurrences  of  several  large  events 
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(M  > 7 with  numerous  aftershocks)  inside  China,  and  extracts  dispersion  infor- 
mation for  a detailed  analysis  on  the  crustal  and  upper  mantle  structures. 

Tectonic  Setting  of  China 

The  tectonic  setting  of  the  four  subplates  of  China  proposed  by  Sun  and 
Teng  (1977)  are: 

1)  The  Chinghai-Tibet  subplate  consists  of  the  area  north  of  the  Himalayas 
including  the  Kun  Lun  Mountains.  Its  eastern  boundary  coincides  with  the 
western  border  of  Szechwan  and  Yunnan  provinces  at  approximately  longitude 

1 05° E.  Its  western  boundary  is  less  well  defined  and  coincides  approximately 
with  the  western  border  of  China.  These  borders  of  the  Chinghai-Tibet  sub- 
plate are  coincident  with  zones  of  high  seismic  activity  and  fossil  sub- 
duction  (Molnar  and  Tapponier,  1975,  1977). 

2)  The  North  China  subplate  (Figure  2),  is  bounded  to  the  north  by  the  Siberian 
platform  (Hamilton,  1970),  to  the  east  by  the  Pacific  Ocean,  to  the  south 

by  the  Tsinling  Shan,  and  to  the  west  by  the  Shansi  graben  and  Ordos  desert. 

The  Tsinling  Shan  consists  of  a series  of  faults  striking  northwesterly  into 
the  Nanshan  fold  belt  and  is  considered  as  both  a distinct  paleontogical  and 
physical  boundary  between  the  two  subplates.  The  Shansi  graben  consists  of  a 
series  of  northeast  striking  faults  which  extend  north  into  the  Ordos  platform. 
These  features  are  evident  by  LANDSAT-1  imagery  as  discussed  by  Cardwell  and 
Isacks  (1976).  A major  part  of  the  North  China  su'oplate  is  the  North  China 
plain,  a sediment-filled  basin  with  an  average  elevation  of  less  than  100 
meters,  surrounded  by  uplifted  areas  of  more  than  1 km  in  elevation.  Although 
the  North  China  subplate  is  considered  an  "intraplate"  area,  it  has  a long 
history  of  seismic  activity  as  indicated  by  the  catalog  of  historical  earth- 
quakes in  China  (Lee,  et  al_.  1976). 

3)  The  South  China  subplate  is  south  of  the  Tsinling  Shan  and  east  of  the 
Tibetan  and  Tsaidam  plateaus.  It  consists  of  the  Yangtze  fault  block  in  the 
northwest  and  the  South  China  fold  block  in  the  southeast  (Sun  and  Teng,  1977). 


Its  eastern  boundary  runs  parallel  to  the  shelf  margin  of  the  southeast 
coast  of  China,  partially  coincident  with  the  subduction  of  the  Pacific 
plate  below  the  Eurasian  plate.  The  presence  of  Melange  in  the  East  Coast 
Mountain  Range  of  Taiwan,  glaucophane  in  the  Central  Mountain  Range  of  Taiwan 
(Sun  and  Teng,  1977),  and  considerable  seismic  activity  in  northeastern  Taiwan 
(Katsumata  and  Sykes,  1969;  Tsai,  et_  al^  1978)  have  been  used  as  evidence  that 
this  region  is  an  active  zone  of  subduction. 

4)  The  Northwestern  China  subplate  consists  of  the  Tarim  and  Turfan  basins, 
the  Tien  Shan  fold  belt,  and  the  area  west  of  the  Ordos  desert.  The  regional 
tectonics  of  this  area  are  not  well  known. 

Data 

The  earthquakes  used  in  the  present  study  (Table  1)  were  chosen  such  that  a 
major  portion  of  their  wave  paths  fall  within  the  four  major  subplates  of  China 
(Fig.  3).  A number  of  earthquakes  are  used  for  each  path  in  order  to  obtain  an 
estimate  on  the  data  repeatability  and  observational  uncertainties.  The  path 
lengths  range  from  approximately  1700  km  to  over  5000  km.  For  the  path  Tangshan 
to  Taipei,  the  accuracy  of  surface  wave  data  for  periods  greater  than  64  seconds 
was  limited.  All  events  were  recorded  at  two  Seismological  Research  Observatory 
(SR0)  stations  located  in  Mashad,  Iran  (MAI0)  and  Taipei,  Taiwan  (TATO). 

SR0  data  in  digital  form  on  tape  is  convenient  for  computer  processing.  The 
SR0  stations  operate  at  high  sensitivity  levels,  thus  expanding  the  detection 
threshold  for  seismic  events  and  making  small -magnitude  events  (M  n,  5 useful  for 
long-  and  short-  period  surface  wave  dispersion  analysis.  Each  seismic  event  must 
be  decoded  and  plotted  from  the  data  tapes.  Pure  Love  and  Rayleigh  waves  are  ob- 
tained by  rotating  the  data  coordinates.  An  example  of  a plotted  seismic  event, 
unrotated  and  rotated,  is  given  in  Figs.  4 and  5.  A complete  discussion  of  handling 
SR0  data  can  be  found  in  Peterson,  et  aU  (1976). 

Method 


A multiple-filtering  analysis  (Dziewonski  and  Hales,  1970;  Hermann,  1973; 
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Tung,  1974;  Seekins  and  Teng,  1976)  has  been  applied  to  the  surface-wave  train 
to  obtain  group  velocity  dispersion  curves.  A typical  contoured  result  over 
the  path  Tangshan  to  Taipei  is  shown  for  both  the  Rayleigh  and  the  Love  com- 
ponents in  Figs.  6,  7 and  8.  In  these  figures  the  dotted  line  represents  the 
average  group  velocities  obtained  from  the  multiple  filtering  analysis.  At 
the  low-period  end  (below  20  seconds)  ambiguity  often  arises  in  drawing  contours 
due  to  interference  by  S waves.  The  countour  lines  will  give  several  local  maxima 
and  will  not  allow  a clearcut  determination  of  group  velocity  values.  The  average 


dispersion  curves  for  Rayleigh  and  Love  components  are  plotted  for  all  events  of 
each  path  in  Figs.  9-16.  Higher-mode  data  are  included  when  they  are  observable 
from  the  multifilter  analysis.  The  statistical  properties  over  each  path  are 
listed  in  Tables  2,  3,  4 and  5. 

The  inversion  of  the  surface  wave  dispersion  data  for  determining  earth 
structure  is  based  on  the  Backus-Gilbert  inversion  method  (1967,  1968,  1970). 
Detailed  descriptions  of  this  method  have  been  presented  (Jackson,  1972;  Wiggins, 
1972;  and  Crosson,  1976).  The  basic  logic  of  the  inversion  theory  and  some  ex- 
perience in  its  numerical  applications  will  be  briefly  discussed  as  follows. 

The  observations,  Y^,  are  related  to  certain  model  parameters,  Xj,  in  some 
known  way,  = A(Xi,X2,  . . . Xp).  A quasi-linear  relationship  is  established 
by  taking  a Taylor  series  expansion  about  some  initial  value  X"0- 
. 3A. 

Y.  = A.(Xo)  + — Ky  • I AX  + higher  order  terms  (1) 

1 1 3Xj  Xo  j 

where  Y . = A^  (X0)  + AY^ 

Ignoring  the  higher  order  terms,  we  obtain  the  functional  relationship  in  matrix 
form: 

AAX  = AY  (2) 

where  AY^  = observed  data  minus  calculated  observations,  and  A is  a matrix  whose 

elements  are  the  partial  derivatives  in  Equation  1.  AX  is  a p x 1 vector  whose 

> 

elements  are  corrections  to  the  initial  model  X0.  AY  is  the  difference  between 
the  observed  and  theoretical  data. 
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In  the  present  case  the  observations  are  group  velocities  of  Love  and  Rayleigh 
waves,  and  the  model  parameters  are  shear  velocities  and  densities  for  a horizon- 
tally layered  earth.  The  partial  derivatives  in  equation  1 are  generated  numeri- 
cally by  the  method  described  by  Rodi  et  al_.  (1975). 

The  classic  least  squares  solution  to  equation  2 is  by  minimizing  the 
Euclidean  length  of  AAX-AY,  or 

II  AAX-AY  II  = (AAX-AY)T(AAX-AY)=  e 

Given  equation  2,  the  solution  is  (Hanson  and  Lawson,  1975): 

A AX  = AY 

ataax  = atay 

AX  = (ATApATAY  . . . 

A 

where  AX  is  the  estimated  correction  vector.  '*“• 

A suitable  inverse  to  this  problem  can  be  obtained  by  using  the  singular 
value  decomposition  (SVD)  of  the  matrix  A. 

A = UAVT 

where  the  columns  of  U are  the  eigenvectors  associated  with  the  columns  of  A, 
the  rows  of  V are  the  eigenvectors  associated  with  the  rows  of  A,  and  A is  a 
diagonal  matrix  of  non-zero  eigenvalues  of  the  matrix  A. 

\\  • * • • 

0 X2 . • • 

n 


0 . . . X, 
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If  we  let 

H = (AW  (3a) 

then 

HAAX  = HAY 

and 

AX  = HA  X 

and 

AX  = HAY 

Using  the  SVO  we  obtain  (Jackson,  1972): 

AX  = (AV'VjAY  (3d) 

H = AV_V 

Here,  H is  known  as  the  pseudo-inverse  of  equation  2. 

The  matrix  product  HA=R  is  known  as  the  resolution  matrix,  and  is  a measure 
of  the  uniqueness  of  the  solution  (Jackson,  1972).  The  rows  of  R are  called  the 

/V  A 

resolving  kernels.  From  equation  3c  it  can  be  seen  that  an  element  of  AX,  AX^, 
may  be  interpreted  as  convolving  the  i ^ row  of  R with  the  vector  AX.  Therefore, 
AX^  can  be  considered  the  weighted  sum  of  nearby  values  (Jackson,  1972). 

Because  ATA  is  nearly  singular,  a number  of  problems  arise:  (1)  The  solution 
vector  becomes  large,  which  can  cause  the  problem  to  leave  the  region  of  linearity; 
(2)  the  solution  oscillates  with  each  iteration.  By  looking  at  the  pseudo-inverse 
H we  can  see  the  cause  of  this  instability. 


(3b) 


(3c) 


and  the  - * s are  the  eigenvalues  of  the  matrix  ATA.  If  A^A  is  near  singular, 
cr.e  or  more  of  the  eigenvalues  will  be  approaching  zero.  From  equation  3 it  can 
be  seen  that  a small  eigenvalue  will  cause  a large  change  in  one  or  more  values 

A 

of  the  correction  vector  AX.  Similarly,  it  has  been  shown  that  the  variance  of 
the  model  parameters  are  inversely  proportional  to  X (Jackson,  1972).  Thus,  small 
X's  produce  large  standard  deviations. 

There  are  two  main  approaches  to  stabilizing  the  inversion  process.  One 


method  is  to  examine  the  eigenvalue  spectrum  of  the  matrix  A.  Small  eigenvalues 
are  removed  when  the  variance  becomes  too  large,  but  removing  the  eigenvalues 
degrades  the  resolution.  This  trade-off  between  resolution  and  variance  has  been 
discussed  by  a number  of  authors  (Jackson,  1972;  Wiggins,  1972;  and  Crosson,  1976). 
Braille  and  Keller  (1975)  used  this  method  for  the  inversion  of  group  velocity 
data. 

Inversion  stability  can  also  be  achieved  by  using  Harquardt's  method 
(Marquardt,  1963)  or  the  stochastic  inverse  (Franklin,  1970).  These  procedures 
lead  to  the  suppression  of  small  eigenvalues.  Franklin's  stochastic  inverse 
estimates  the  parameter  in  the  presence  of  noise. 


Where  n is  a vector  of  observational  noise,  Franklin  (1970)  shows  that  the 
solution  to  equation  4 is  given  by: 

AX  = WAT(AWAT  + ED)" *AY  (5) 

where  W is  the  covariance  matrix  of  the  parameter  and  ED  is  the  covariance  matrix 
of  the  observations. 

Marquardt's  method  considers  a minimization  of  the  functional 
(AAX  - AY)T(EDTED)  (AAX  - AY)  - AXTIWTW)  AX  (6) 

with  the  solution  given  by: 

AX  = (ATED_1  A + W "VVed"1  y (7) 

T - 1 T 

Letting  A ED  = A , equation  7 becomes: 

AX  = (ATA  + a2I)  ATAY  (8) 

2 

where  a is  the  variance  of  the  parameter. 

By  looking  at  the  SVD  of  equation  8 we  can  see  how  Marquardt's  method  tapers 
the  eigenvalue  spectrum: 

AX  = V|  (A2  + a2I)-1A|UTAY  (9) 

An  element  in  brackets  from  equation  9 equals:  (Crosson,  1976). 

I (A  + a I)"  A|.  = \./[\  ■ + o ) 

So  as  A.j  approaches  zero,  the  parameter  AX^  also  approaches  zero.  A similar 
analysis  for  the  stochastic  inverse  shows  that  it  too  suppresses  small  eigenvalues. 


The  quantity  a can  also  be  thought  of  as  a tradeoff  parameter  between 
resolution  and  variance  (Der  «rt  al_-»  1970;  Wiggins,  1972;  Crosson,  1976).  In 
Marquardt's  method,  the  resolution  equals  (Crosson,  1976): 

R = HA  = (ATA  + o2I)ATA  = V|(A  + X2I)-1A2 | VT  (10) 


For  the  stochastic  inverse 


2 2 1 2 t 

R * | (A  + o I)'  A VV1 


The  quantity  in  brackets  of  equations  10  and  11  is  a diagonal  matrix  whose 
elements  have  the  form: 


2 2 2 

A.  / X + a ) 

Z 2 

If  z =0,  then  R - I,  but  as  a approaches  zero,  then  the  variance  of  the 

2 

parameters  becomes  larger.  We  then  adjust  the  size  of  o until  there  is  an 
acceptable  trade-off  between  standard  deviation  and  resolution. 

Two  other  elements  that  must  be  analyzed  in  the  inversion  process  are  the 
apriori  (SDX0)  and  postpriori  (SDX)  standard  deviations  of  the  inversion  para- 
meters: 

SDX0  = DT(oWP"1)D 

SDX  = v^(ATED-1A  + oWP-1pD 

where  D = a row  vector  of  a delta  matrix  = (0,0,. . . ,1 ,0,. . . ,0) . This  analysis  will 
give  a measure  of  how  much  new  information  is  coming  from  the  data.  If  SDX0  and 
SDX  are  approximately  the  same,  we  know  that  very  little  additional  information 
is  coming  from  the  data. 

The  advantages  of  Marquardt's  method  and  the  stochastic  inverse  are  that 
(1)  the  SVD  of  A is  not  explicitly  determined,  and  (2)  there  is  no  decision  to  be 
made  on  the  rank  of  the  matrix  A^A. 
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I A combination  of  these  two  methods  is  used  in  our  inversion  process.  When 

the  number  of  parameters  is  greater  than  the  number  of  data,  the  stochastic 
inverse  is  used.  If  the  number  of  observations  is  greater  than  the  number  of 
parameters,  Marquardt's  method  is  used.  This  procedure  leads  to  the  minimum 
number  of  computations. 

i ; 

Results  and  Discussion 

Initial  models  for  the  inversion  were  based  on  the  results  of  a previous 
surface  wave  study  done  by  Tung  (1974).  Starting  models  for  all  paths  were 

I [ 

initially  inverted  for  both  shear  velocity  and  density.  During  the  inversion, 

j | 

large  instabilities  of  the  density  parameter  were  found  in  the  crust.  This  was 

, 

| corrected  by  changing  the  starting  model  and  inverting  for  the  compressional 

wave  velocity  with  the  shear  velocity  and  density.  As  discussed  by  Wiggens  (1972), 
Der,  et  a£.  (1978),  and  Jackson  (1976),  the  partial  derivatives  for  shear  velocity 
are  larger  than  those  for  compressional  velocity  or  density.  Therefore  the  reso- 
lution for  shear  wave  velocity  is  considerably  better  than  the  other  two  para- 
meters. Only  the  resolution  matrix  for  the  partial  derivatives  with  respect  to 
shear  velocity  shall  be  presented. 

1.  Tibetan  Plateau 

Two  separate  paths  across  the  Tibetan  plateau  were  used  (Fig.  3).  Both 
paths  originate  from  earthquakes  located  on  the  eastern  boundary  of  the  Tibetan 
plateau.  The  first  group  occurred  in  Szechwan  province  and  the  second  group 
occurred  to  the  south  in  Yunnan  Province.  A list  of  the  surface  waves  is  presented 
in  Table  6.  For  each  path  across  the  Tibetan  plateau  inversions  were  performed 
using  the  fundamental  and  first  higher  mode  Rayleigh  waves.  A fit  of  the 
theoretical  model  to  the  observed  group  velocity  is  presented  in  Fig.  17  and  18. 

A comparison  between  the  average  continental  dispersion  for  Rayleigh  waves  (Ewing, 
et  £^,  1957)  and  Rayleigh  waves  for  the  four  paths  across  China  are  presented  in 
Fig.  19.  The  group  velocities  across  the  Tibetan  plateau  are  substantially  lower 


than  the  group  velocities  for  the  average  continental  structure,  with  differences 
extending  to  periods  of  120  seconds,  indicating  abnormal  upper  mantle  structure. 

Group  velocities  across  the  Tibetan  plateau  are  the  lowest  of  all  the 
regions  of  China.  Rayleigh  wave  group  velocities  ranged  between  2.60  to  3.65  km/sec 
and  steep  slopes  occurred  in  the  dispersion  curves  between  40  to  60  seconds,  where 
the  change  in  group  velocity  was  as  high  as  0.6  km/sec.  Love  waves  were  not 
used  in  the  inversions  due  to  the  absence  of  amplitude  maxima  in  the  dispersion 
curves  between  30  to  60  seconds.  Earlier  studies  (Thatcher  and  Brune,  1969;  James, 
1971)  have  discussed  the  problem  of  interference  of  Love  waves  crossing  inhomogeneous 
media.  Tung  (1974)  attributed  this  problem  on  the  Tibetan  plateau  to  Love  waves 
following  a non-least  time  path.  A possible  structural  implication  is  that  there 
may  be  structures  along  the  path  corresponding  to  wave  lengths  of  the  Love  waves, 
causing  the  wave  energy  to  travel  along  a non-least  time  path  with  slower  velocities. 

Starting  and  final  shear  wave  velocity  models  with  standard  deviation  bars 


for  the  path  Szechwan  to  Mashad  are  presented  in  Figure  20.  Table  7 presents 
the  initial  and  final  model  parameters  and  the  resolution  matrix  for  each  model. 
The  resolution  matrix  for  the  eight-layer  model  is  presented  in  Table  8.  From 
the  width  of  the  resolving  kernels  it  is  possible  to  resolve  layers  of  15  km  in 
the  upper  crust,  15-20  km  in  the  middle  crust,  and  25-30  km  in  the  lower  crust. 

Starting  and  final  shear  wave  velocity  models  with  standard  deviation  bars 
for  the  path  Yunnan  to  Mashad  are  presented  in  Figure  21.  Tables  9-10  present 
the  initial  and  final  model  parameters  and  the  resolution  matrix.  From  the  width 
of  the  resolving  kernels  it  is  possible  to  resolve  a 10-15  km  layer  at  the  top 
of  the  crust,  15-20  km  layer  in  the  middle  of  the  crust,  and  a 20-25  km  layer 
at  the  base  of  the  crust. 

Resolution  in  the  mantle  is  poor  for  both  paths.  The  standard  deviations 
were  found  to  be  less  than  0.15  km/sec  for  all  the  models  presented.  Increasing 
the  resolution  in  the  crust  and  mantle  leads  to  both  large  standard  deviations 
and  instabilities  in  the  shear  wave  velocity  model  (Rosenthal,  1977).  A measure 
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of  the  amount  of  new  information  obtained  from  the  data  is  given  by  comparing 
the  magnitudes  of  the  apriori  standard  deviation  (SDX0)  to  the  postpriori  standard 
deviations  (SDX)  given  in  Tables  8 and  10.  In  the  upper  mantle  there  is  only  a 
small  amount  of  new  information  being  obtained  from  the  data,  and  the  final 
model  is  most  affected  by  the  starting  model.  As  a result,  a low  velocity  layer 
in  the  upper  mantle  can  only  be  inferred  from  the  small  amount  of  new  information 
being  obtained.  Additional  data  for  periods  longer  than  120  seconds  and  good  Love 
wave  data  would  significantly  increase  the  amount  of  new  information. 

The  shear  wave  velocity  structure  for  both  paths  crossing  the  Tibetan  plateau 
is  similar.  The  crust  is  approximately  70  km  thick  and  may  be  broken  into  four 
layers.  The  first  layer  has  shear  velocities  that  range  between  2.87-3.06  km/sec, 
suggesting  a layer  approximately  15  km  thick  of  partly  sedimentary  origin.  Tung 
(1974)  and  Mu,  et^  aJL  (1974)  discussed  the  existence  of  thick  sedimentary  layers 
in  the  Mount  Jolmo  Lungma  (Mount  Everest)  region  of  the  Himalayas.  Both  Birch 
(1963)  and  Simmons  (1964)  reported  shear  velocities  for  sedimentary  rocks  at  10 
to  15  km  depth  of  approximately  3.00  to  3.08  km/sec  A second  layer  has  shear 
wave  velocities  ranging  from  3.20  to  3.40  km/sec  end  is  approximately  15  km  thick, 
suggesting  a layer  of  granitic  or  equivalent  composition.  A third  layer  has 
shear  velocities  between  3.40  to  3.70  km/sec,  is  approximately  15  to  20  km  thick, 
and  also  seems  to  consist  of  granite.  For  comparison,  Birch  and  Bancroft  (1958), 
Birch  (1963)  and  Simmons  (1964)  reported  shear  velocities  for  granite  under  3 to 
5 kilobars  of  pressure  of  3.45-3.70  km/sec.  Simmons  (1964)  reported  that  granite 
with  a shear  wave  velocity  of  3.56  km/sec  at  4 kilobars  pressure  has  a shear  wave 
velocity  of  3.79  km/sec  at  pressures  of  10  kilobars.  The  fourth  layer  with  shear 
wave  velocities  between  3.85  to  3.95  km/sec  corresponds  to  granite  at  pressures 
of  approximately  20  kilobars,  which  would  correspond  to  depths  of  60  to  70  km. 

The  shear  wave  velocity  profile  for  the  Tibetan  plateau  is  in  agreement 
with  those  values  presented  by  Tung  (1974)  and  Bird  and  Toksoz  (1977).  Chun 
and  Yoshi  (1977),  by  inverting  Rayleigh  and  Love  waves  for  many  different  paths 
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across  the  Tibetan  plateau,  reported  a crustal  thickness  of  approximately  70  km 
and  observed  a velocity  reversal  in  the  middle  crust.  Due  to  a 15  to  20  km 
resolution  in  the  middle  crust,  it  is  doubtful  a velocity  reversal  could  be 
observed  for  data  over  the  two  paths  presented. 

Due  to  the  small  amount  of  information  (apriori  vs.  postpriori  standard 
deviations)  derived  from  inverting  Rayleigh  waves  in  the  upper  mantle,  the 
above  structural  conclusions  are  only  tentative.  The  upper  mantle  shear  velocities 
are  lower  than  either  the  Gutenburg  earth  model  of  Tacheuchi,  et  aK  (1964)  or 
the  Brune  and  Dorman  (1963)  model  for  the  Canadian  Shield  (Fig.  22).  The  velocity 
at  the  top  of  the  mantle  is  approximately  4.4  km/sec,  with  indications  of  a low 
velocity  layer  at  90  km  depth.  Since  the  low  velocity  layer  can  only  be  inferred, 
the  existence  of  extremely  low  velocities  of  4.10  km/sec  is  tentative.  Low 
shear  velocities  in  the  upper  mantle  of  Tibet  have  been  reported  by  Tung  (1974), 
which  suggests  a high  degree  of  partial  melting  and  a low  resistance  to  movement 
(Ancerson,  et  al_.  1972). 

2.  North  and  Northwestern  China 

The  path  between  Tangshan  and  Mashad,  Iran  (Fig.  3)  crosses  the  North  China 
subplate  and  the  Northwestern  China  subplate  including  the  Tien  Shan  fold  belt. 

Tung  (1974)  found  crustal  thicknesses  of  30,  40  and  45  km  for  Northeast,  Central 
and  Northwest  China  respectively.  Crustal  thicknesses  in  the  Tien  Shan  and  Pamir 
Mountains  have  been  reported  as  high  as  65  km  (Arkhangel 'skaya,  1964). 

A list  of  the  earthquakes  and  surface  wave  trains  used  over  the  path  Tangshan 
to  Mashad,  Iran  are  listed  in  Table  11.  The  fit  of  the  theoretical  to  observed 
group  velocities  for  the  Rayleigh  and  Love  waves  is  given  in  Fig.  23.  The 
theoretical  Rayleigh  wave  group  velocity  is  lower  than  the  average  continental 
Rayleigh  wave  dispersion  (Ewing,  et  aK  1957).  Starting  and  final  shear  wave 
velocity  models  with  standard  deviation  bars  are  presented  in  Fig.  24.  Tables 
12  and  13  present  the  starting  and  final  model  parameters  and  the  resolution 
matrix.  From  the  width  of  the  resolving  kernels  presented  in  Table  13,  it  is 
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When  Rayleigh  and  Love  waves  are  inverted  simultaneously,  a considerable 
amount  of  new  information  may  be  derived  from  the  data.  Observing  the  magni- 
tude of  the  apriori  (SDX0)  vs.  the  postpriori  (SOX)  standard  deviations,  a con- 
siderable amount  of  new  information  is  obtained  at  depths  of  100  km  (seventh 
layer).  Below  this  depth,  less  new  information  is  obtained  from  the  data,  and 
the  final  model  is  most  affected  by  the  starting  model. 

The  average  crustal  thickness  over  the  path  from  Tangshan  to  Mashad  is 
approximately  45  km.  Since  this  path  is  not  considered  a pure  path,  only  average 
shear  velocities  are  presented  (Fig.  3).  A three-layer  crust  seems  to  be  an 
adequate  model  for  this  path.  The  first  layer  is  approximately  10  km  thick  and 
has  an  average  shear  velocity  of  approximately  3.00  km/sec,  suggesting  a layer 
at  least  partly  sedimentary  (Lirch,  1963;  and  Simmons,  1964).  The  second  layer 
is  15  to  20  km  thick  and  has  an  average  shear  wave  velocity  of  3.55  km/sec.  The 
layer  at  the  base  of  the  crust,  20  to  25  km  thick,  has  an  average  shear  velocity 
of  approximately  3.80  km/sec,  suggesting  a layer  of  granitic  origin  (Birch,  1963; 
Simmons,  1964).  The  average  shear  wave  velocity  at  the  top  of  the  mantle  is  4.40 
to  4.45  km/sec  and  the  low  velocity  layer  begins  at  approximately  90  km.  Using  a 
three- layer  earth  model  proposed  by  Dorman  (1959),  Schechkov  (1961)  was  able  to 
fit  mixed  path  dispersion  data  for  Northwest  China  to  a theoretical  three-layer 
model  of  average  crustal  thickness  of  50  km:  the  first  layer  has  a shear  velocity 
of  3.40  km/sec  and  is  20  km  thick,  the  second  layer  has  a 3.53  km/sec  shear  velocity 
and  is  30  km  thick,  and  a semi-infinite  third  layer  has  a 4.50  km/sec  shear  velocity 
for  the  mantle.  Applying  the  concept  of  regionalization  to  paths  crossing  North- 
west China,  Tung  (1974),  presented  a three-layer  crust  of  45  km  thickness.  Shear 
velocities  calculated  by  Tung  (1974)  for  the  crust  and  upper  mantle  are  lower  than 
the  mixed  path  average  shear  velocities  presented.  The  low  velocity  zone  presented 
is  not  as  prominent  as  Tung  (1974)  reported,  and  can  only  be  inferred  from  the 


Information  obtained  from  the  data. 

A comparison  of  the  Canadian  Shield  model  (Tacheuchi,  et  aK  1964)  with 
the  average  shear  velocity  model  for  Tangshan  to  Mashad  is  presented  in  Fig. 

25.  Both  models  have  higher  shear  velocities  in  the  crust  and  at  the  top  of 
the  mantle  than  the  theoretical  model  presented  for  the  path  Tangshan-Mashad. 

The  tectonics  over  the  path  between  Tangshan-Mashad  may  account  for  the  lower 
velocity  at  the  top  of  the  mantle.  The  catalog  of  historical  earthquakes  pre- 
sented by  the  Chinese  Institute  of  Geophysics  Academia  Sinica  (Lee,  et^  al_.  1976) 
shows  considerable  seismic  activity  throughout  the  regions  traversed  by  the 
surface  waves.  Successive  stages  of  subduction  of  the  Indian  plate  below  the 
Eurasian  plate  are  believed  to  have  taken  place  in  areas  as  far  north  as  the 


Tarim  basin  and  Mongolia  (Chang  and  Zeng,  1973;  Zonenshain,  1973).  These  are 
consistent  with  a thicker  crust,  lower  group  velocities  and  lower  shear  velocities 
in  the  upper  mantle.  These  results  suggest  that  the  origin  of  the  Tien  Shan  fold 
belt  may  also  be  a result  of  the  collision  of  the  Eurasion  and  Indian  plates. 

3.  Coastal  China 


For  most  of  the  path  between  Tangshan  to  Taipei  the  surface  waves  are 
parallel  to  the  east  coast  of  China.  According  to  the  tectonic  map  of  the 
Institute  of  Geology  of  The  Chinese  Academy  of  Sciences  (from  which  Fig.  1 is 
derived),  the  crustal  thickness  along  most  of  the  path  is  about  30  km.  A starting 
model  for  inversion  was  chosen  to  reflect  the  30  km  crustal  thickness.  A list 


of  the  earthquakes  and  surface  wave  trains  used  in  the  analysis  is  given  in  Table 
14.  The  fit  of  the  theoretical  to  the  observed  Rayleigh  and  Love  wave  group 
velocities  is  presented  in  Fig.  26.  A comparison  between  the  theoretical  Rayleigh 
and  the  average  continental  Rayleigh  wave  group  velocities  (Ewing,  et  aK  1957) 
shows  lower  group  velocities  of  both  short  and  intermediate  periods. 

The  starting  and  final  shear  wave  velocity  models  with  standard  deviation 
bars  are  presented  in  Fig.  7.  The  shear  velocity  parameters  and  resolution  matrix 
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are  presented  in  Tables  15  and  16.  From  the  width  of  the  kernels  of  the  resolution 
matrix  (Table  16),  it  is  possible  to  resolve  a 5 to  7 km  layer  at  the  top  of  the 
crust,  a 10  km  layer  in  the  middle  of  the  crust,  and  a 15  to  20  km  layer  at  the 
base  of  the  crust.  An  examination  of  the  apriori  (SDX0)  and  postpriori  (SDX) 
standard  deviations  indicates  that  a small  amount  of  information  exists  in  the 
data  below  the  fifth  layer.  Thus  longer  period  data  is  needed  for  increased  in- 
formation in  the  mantle.  Since  the  path  length  from  Tangshan  to  Taipei  is  approxi- 
mately 1700  km,  periods  greater  than  64  seconds  will  not  have  proper  mode  separation. 
Consequently,  it  is  unlikely  that  for  this  path  new  information  below  the  top  of 
the  mantle  can  be  obtained  and  shear  velocity  in  the  mantle  is  most  affected  by 
the  starting  model. 

A crust  with  a thickness  of  30  km  and  composed  of  three  layers  is  an  adequate 
model  for  the  structure  over  the  path.  The  first  layer,  approximately  5 to  7 km 
thick,  seems  to  be  a sedimentary  layer  with  shear  velocity  of  approximately  2.95 
km/sec.  The  second  layer  is  n,  10  km  thick.  The  third  layer,  approximately  15  to 
20  km  thick,  corresponds  to  a shear  velocity  of  approximately  3.85  km/sec.  in  the 
mantle,  a low  velocity  layer  can  only  be  inferred  due  to  the  lack  of  new  information 
from  the  data.  A comparison  of  the  Canadian  Shield  model  (Brune  and  Dorman, 

1953)  and  the  Gutenberg  earth-model  (Tacheuchi,  e^t  al_.  1964)  with  the  shear  velocity 
profile  shows  higher  velocities  in  the  mantle  and  the  base  of  the  crust  (Fig.  28). 

The  shear  velocity  of  4.46  km/sec.  at  the  top  of  the  mantle  is  lower  than  the 
shear  velocity  of  the  Canadian  Shield  (4.65  km/sec.).  Tung  (1974)  also  calculated 
a mantle  shear  velocity  of  4.42  km/sec.  at  the  top  of  the  mantle  for  a path  crossing 
Northeast  China  recorded  in  Hong  Kong.  The  tectonics  over  the  path  from  Tangshan 
to  Taipei  explains  the  lower  velocities  in  the  mantle. 

The  region  traversed  is  considered  an  "intraplate"  region,  although  a con- 
siderable amount  of  seismic  history  is  evident  from  the  catalog  of  historical  earth- 
quakes compiled  by  the  Chinese  Institute  of  Geology  Academia  Sinica  (Lee,  et  at.  1976). 
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For  example,  in  the  northeastern  portion  of  the  South  China  subplate  is  the  Tanlu 
fault,  which  may  be  traced  from  Landsat  photographs  (Lee,  e_t  aK  1976)  northward 
across  the  border  with  the  North  China  subplate,  through  Shantung  Province,  the 
Gulf  of  Chihli  (Pohai)  and  into  Liaoning  Province  (north  of  Tangshan).  The  Tangshan 
earthquake  (July  1976)  and  the  Haichang  earthquake  (February  1975)  indicate  that  the 
northern  portion  of  the  path  from  Tangshan  to  Taipei  passes  through  a tectonically 
active  region. 

The  southern  portion  of  the  path  traverses  a region  affected  by  the  Philippine 
Sea  plate  that  is  subducted  beneath  the  Eurasion  plate.  The  tectonics  of  this 
region  have  been  discussed  by  Sun  and  Teng  (1977).  Considerable  seismic  activity 
is  found  throughout  northern  and  eastern  Taiwan,  particularly  along  a major  north- 
south  trending  fault  paralleling  the  east  coast  of  Taiwan.  Sun  and  Teng  (1977) 
discussed  the  importance  of  this  fault  to  a subduction  zone  of  Mesozoic  or  Cenozoic 
age.  The  eastern  provinces  of  Fuchien,  Kiangsu,  Chiekiang,  Kwantung,  Anhwai  and 
Kiangsi  are  covered  by  large-scale  volcanic  intrusions  which  are  of  Mesozoic  age. 

These  volcanic  intrusions  are  attributed  to  a magma  source  which  was  a consequence 
of  the  melting  from  the  subduction  of  the  Pacific  Sea  plate.  Finally,  higher  than 
normal  heat  flow  (3  to  4 heat  flow  units)  has  been  observed  north  of  Taiwan  which 
is  associated  with  areas  of  plate  subduction  (Yansuii,  et^  al_.  1970).  All  these 
factors  would  contribute  to  a shear  wave  velocity  profile  lower  than  normal  in  the 
upper  mantle. 

Conclusions 

By  using  a multiple-filter  technique,  group  velocities  for  various  regions  of 
China  were  obtained  from  earthquakes  recorded  at  SRO  stations  in  Mashad,  Iran  and 
Taipei,  Taiwan.  For  each  sampling  path,  a number  of  events  were  used  with  practically 
the  same  source-station  geometry  to  derive  the  dispersion  data.  An  estimate  of  un- 
certainty of  the  data  was  thus  obtained.  A nonlinear  least  square  technique  was 
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used  to  invert  the  surface  wave  group  velocity  for  shear  wave  velocity  structures. 
The  findings  concerning  the  nature  of  the  crust  and  mantle  are: 

1.  Group  velocities  obtained  for  the  pure  path  across  the  Tibetan  plateau 
are  unusually  low  when  compared  with  the  average  continental  dispersion.  Results 
of  the  inversion  indicate  that  a 4-layer  70-km  thick  crust  is  an  adequate  model 
for  the  Tibetan  plateau. 

2.  Group  velocities  for  the  mixed  path  between  Tangshan  and  Mashad,  Iran 
are  lower  than  the  average  continental  dispersion.  Results  of  the  inversion 
indicate  that  a 3-layer,  45  km  thick  crust  is  an  adequate  model. 

3.  Group  velocities  for  the  path  between  Tangshan  and  Taipei,  Taiwan  are 
closest  of  the  three  regions  to  average  continental  dispersion.  Results  of  the 
inversion  indicate  that  a 3-layer  30-km  crust  is  an  adequate  model. 

For  all  three  models,  the  structure  of  the  upper  mantle  from  information 
provided  by  the  data  is  tentative.  In  comparing  the  apriori  (SDX0)  to  the  post- 
priori  (SDX)  standard  deviations  in  the  mantle,  it  was  observed  that  the  final 
model  was  influenced  by  the  starting  model  more  than  by  the  data.  Attempts  to 
increase  the  information  from  the  data  and  the  resolution  result  in  either  un- 
acceptably large  standard  deviations  or  unstable  shear  velocity  models.  Only  in 
the  case  of  data  for  the  path  from  Tangshan  to  Mashad,  for  which  inversions  were 
performed  simultaneously  for  long  period  Rayleigh  and  Love  waves,  was  new  infor- 
mation obtained  about  the  upper  mantle. 

Previous  surface  wave  studies  of  China  have  applied  a trial -and-error  inversion 
approach  to  obtain  a shear  velocity  model  from  Rayleigh  wave  group  velocity  data. 
Proper  consideration  has  not  been  given  to  either  the  information  obtained  from 
the  data,  or  the  resolution.  In  all  these  studies  the  influence  of  the  starting 
model  upon  the  final  model  was  not  discussed,  thus  the  validity  of  the  estimated 
upper  mantle  shear  velocities  presented  is  rather  questionable. 

Both  Rayleigh  and  Love  waves  should  be  inverted  simultaneously  to  obtain 
more  information  from  the  data.  This  may  not  always  be  possible  due  to  crustal 


anisotrophy.  Higher  mode  data  are  very  important,  but  generally  the  ex- 
citation of  higher  modes  is  not  strong  for  the  shallow  earthquakes  used 
in  this  study.  In  this  regard,  future  large  earthquakes,  in  Tibet  for 
instance,  of  intermediate  to  deep  focal  depths  would  be  important  sources 
for  higher  mode  data  and  they  warrant  special  attention.  Further,  the 
dispersion  curves  must  not  be  contaminated  by  surface  waves  traveling 
non-least  time  paths  or  by  complex  structure  along  the  path.  This  often 
presents  a problem  in  tectonically  active  regions,  resulting  in  the 
inability  to  obtain  a good  Love  wave  dispersion  curve  from  the  data. 
Finally,  the  problems  due  to  epicentral  distance  must  be  considered.  If 
the  epicentral  distance  is  too  short,  proper  surface  wave  mode  separation 
will  not  be  achieved.  Only  after  those  factors  are  considered  can  one 
arrive  at  a better  shear  wave  velocity  structure  by  means  of  surface  wave 
inversion  for  a specific  geological  province. 
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TABLE  1 

Earthquakes  and  Recording  Stations 


Used  in  this  Studv 


Location 

39.8°  N 
118.6°  E 
N.E.  China 

40.1°  N 
118.3°  E 
N.E.  China 

39.9°  N 
118.0°  E 
N.E.  China 

39.8°  N 
117.8°  E 
N.E.  China 

39.5°  N 
117.9°  E 
N.E.  China 

39.7°  N 
118.5°  E 
N.E.  China 

40.2°  N 
118.9°  E 
N.E.  China 

39.6°  N 
118.5°  E 
N.E.  China 

32.753°  N 
104.157°  E 
Szechwan 

32.393°  N 
104.139°  E 
Szechwan 

32.571°  N 
104.152°  E 
Szechwan 


Origin 

Time 

July  28,  1976 

15  hr  35  min 

55.3  sec 

July  28,  1976 

16  hr  24  min 

04.7  sec 

July  29,  1976 
1 hr.  01  min 

04.1  sec 

July  30,  1376 

21  hr  23  min 

13.8  sec 

Aug  2,  1976 
9 hr  16  min 
00.5  sec 

Aug  8,  1976 
1 hr  0.9  min 

12.4  sec 

Aug  8,  1976 

22  hr  41  min 
34.3  sec 

Aug  14,  1976 
16  hr  02  min 

44.5  sec 

Aug  16,  1976 
14  hr  6 min 

45.9  sec 

Aug  19,  1976 
12  hr  49  min 
47.7  sec 

Aug  21,  1976 
21  hr  49  min 

54.2  sec 


Recording  Station 


MAIO 


MAIO 


MAIO 


MAIO 


MAIO 


MAIO 


MAIO 


MAIO 


MAIO 


MAIO 


MAIO 


Magnitude 


TABLE  1 (Continued) 


Location 

Origin 

Time 

Recording  Station 

Magnitude 

32.492°  N 
104.181°  E 
Szechwan 

Aug  23,  1976 

3 hr  30  min 

7.6  sec 

MAIO 

6.2 

32.460°  N 
104.152°  E 
Szechwan 

Sept  1,  1976 

1 hr  6 min 

51.8  sec 

MAIO 

5.1 

24.343°  N 
98.642°  E 
Yunnan 

May  31,  1976 

5 hr  8 min 

28.5  sec 

MAIO 

5.5 

24.191°  N 
98.676°  E 
Yunnan 

July  3,  1976 

16  hr  33  min 
23.1  sec 

MAIO 

5.3 

39.3°  E 
113.6°  E 
N.E.  China 

July  28,  1976 
15  hr  35  min 
55.3  sec 

TATO 

5.4 

39.9°  N 
118.8°  E 
N.E.  China 

July  29,  1976 

1 hr  01  min 
04.1  sec 

TATO 

5.1 

39.8°  N 
118.9°  E 
N.E.  China 

July  30,  1976 

TATO 

5.4 

39.6°  N 
117.8°  E 
N.E.  China 

Aug  1,  1976 

20  hr  53  min 
53.6  sec 

TATO 

4.6 

39.9°  N 
118.8°  E 
N.E.  China 

Sept  6,  1976 

17  hr  02  min 
01.5  sec 

TATO 

4.8 

MAIO:  Mashad,  Iran  36.30°  N 

59.49°  E 

TATO:  Tapei,  Taiwan  24.976°  N 

121.489°  E 


■ 


64.0 

3.49 

+ .06 

3.85 

+ .15 

51.2 

3.48 

+ .09 

3.71 

+ .10 

42.7 

3.46 

+ .11 

3.69 

+ .06 

36.6 

3.48 

+ .11 

3.59 

+ .05 

32.0 

3.30 

+ .09 

3.52 

+ .08 

28.4 

3.29 

+ .05 

3.37 

+ .05 

25.6 

3.15 

+ .08 

3.31 

+ .05 

23.3 

3.04 

+ .05 

3.19 

+ .05 

21.3 

3.00 

+ .06 

3.19 

+ .05 

19.7 

2.84 

+ .06 

3.12 

+ .09 

18.3 

2.76 

+ .09 

3.12 

+ .07 

17.1 

2.72 

+ .06 

3.04 

+ .09 

16.0 

2.67 

+ .05 

3.02 

+ .04 

15.1 

2.63 

+ .10 

2.93 

+ .05 

13.5 

2.58 

+ .13 

2.89 

+ .10 

12.2 

2.52 

+ .10 

2.85 

+ .15 

11.1 

2.59 

+ .18 

2.90 

+ .07 

10.2 

2.64 

+ .17 

2.91 

+ .15 

1 " 


TABLE  3 


Tangshan-Mashad  Observed  Dispersion  Data 


I 


■ 


Period (sec) 


Rayleigh  Wave 
(krn/sec) 


Love  Wave 
(km/sec) 


102.4 

3.73 

+ .20 

4.00 

+ .20 

85.3 

3.65 

+ .18 

3.91 

+ .13 

73.1 

3.69 

+ .10 

3.91 

+ .10 

64.0 

3.57 

+ .10 

3.85 

+ .05 

56.9 

3.65 

+ .18 

3.74 

+ .06 

51.2 

3.54 

+ .05 

3.70 

+ .11 

46.5 

3.47 

T .05 

3.66 

+ .10 

42.6 

3.43 

+ .05 

3.56 

+ .06 

39.4 

3.35 

+ .05 

3.46 

+ .05 

36.6 

3.22 

+ .05 

3.38 

+ .05 

34.1 

3.14 

+ .04 

3.36 

+ .05 

32.0 

3.10 

+ .05 

3.31 

+ .05 

30.1 

3.01 

+ .06 

3.29 

+ .06 

26.9 

2.93 

+ .06 

3.  21 

+ .05 

24.4 

2.91 

+ .08 

3.18 

+ .05 

22.3 

2.85 

+ .10 

3.16 

+ .07 

20.5 

2.74 

+ .06 

3.17 

+ .05 

18.9 

2.78 

+ .10 

3.14 

+ .05 

16.0 

2.82 

+ .18 

3.06 

+ .06 

13.8 

2.75 

+ .07 

3.01 

+ .06 

12.2 

2.79 

+ .08 

2.95 

+ .10 

10.9 

2.69 

+ .08 

2.89 

+ .15 

TABLE  4 


Yunnan-Mashad  Observed  Dispersion  Data 


Rayleigh  Wave (km/sec) 

Period  (sec)  Fundamental First  Higher  Mode 


102.4 

3.52 

+ 

.15 

85.3 

3.49 

+ 

.15 

73.1 

3.41 

+ 

.10 

64.0 

3.25 

+ 

.15 

56.9 

3.12 

+ 

.10 

51.2 

3.05 

+ 

.10 

46.5 

3.00 

+ 

.10 

42.7 

2.95 

+ 

.10 

39.4 

2.87 

+ 

.07 

36.6 

2.79 

+ 

.07 

34.1 

2.73 

+ 

.07 

32.0 

2.71 

+ 

.07 

30.1 

2.73 

+ 

.15 

27.0 

2.67 

+ 

.15 

24.4 

2.85 

+ 

.15 

22.3 

2.79 

+ 

.05 

20.5 

2.79 

+ 

.05 

4.09 

+ .15 

18.9 

2.77 

+ 

.05 

3.80 

+ .17 

16.0 

2.77 

+ 

.06 

3.48 

+ .10 

13.8 

2.77 

+ 

.06 

12.2 

2.64 

+ 

.15 

3.17 

+ .08 

10.9 

2.71 

+ 

.10 

3.15 

+ .07 

TABLE  5 


Szechwan-Mashad  Observed  Dispersion  Data 


Rayleigh  Wave (km/sec) 

Period (sec)  Fundamental  First  Higher  Mode 


102.4 

3.54 

+ 

.17 

85.3 

3.56 

+ 

.09 

73.1 

3.49 

+ 

.06 

64.0 

3.41 

+ 

.07 

56.9 

3.28 

+ 

.06 

51.2 

2.96 

+ 

.15 

46.5 

2.96 

+ 

.15 

42.7 

2.95 

+ 

.07 

39.4 

2.87 

+ 

.06 

36.7 

2.87 

+ 

.08 

34.1 

2.87 

+ 

.09 

32.0 

2.79 

+ 

. 10 

30.1 

2.76 

+ 

.09 

26.9 

2.75 

+ 

.06 

24.4 

2.74 

+ 

.07 

22.3 

2.76 

+ 

.04 

20.5 

2.69 

+ 

.05 

18.9 

2.71 

+ 

.04 

16.0 

2.73 

+ 

.04 

3.71 

+ .20 

13.8 

2.70 

+ 

.07 

3.45 

+ .15 

12.2 

2.73 

+ 

.13 

3.35 

+ .15 

10.8 

3.16 

+ .10 

Aug  16, 

1976 

MAIO 

Vertical 

Rayleigh 

4080 

Aug  1 6 , 

1976 

MAIO 

Radial 

Rayleigh 

4080 

Aug  19 , 

1976 

MAIO 

Vertical 

Rayleigh 

4078 

Aug  19 , 

1976 

MAIO 

Radial 

Rayleigh 

4078 

Aug  21, 

1976 

MAIO 

Vertical 

Rayleigh 

4095 

Aug  2 1 , 

1976 

MAIO 

Radial 

Rayleigh 

4095 

Aug  23, 

1976 

MAIO 

Vertical 

Rayleigh 

4088 

Aug  2 3 , 

1976 

MAIO 

Radial 

Rayleigh 

4088 

Sept  1, 

1976 

MAIO 

Vertical 

Rayleigh 

4090 

Sept  1, 

1976 

MAIO 

Radial 

Rayleigh 

4090 

TABLE  8 : Resolution  Matrix  fof  Path  Szechwan-Mashad  (8  Layers) 
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TABLE  9 : Model  Fit  for  Path  Yunnan-Mashad 
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Seismogram  Components  Used  for  the 
Path  Tanqshan  to  Mashad 


Earthquake 

Recording 

Station 

Component 

Wave  Type 

Distance 

(km) 

July 

29. 

1976 

KAIO 

Vertical 

Rayleigh 

5119 

July 

29, 

1976 

MAIO 

North-South 

Love 

5119 

July 

30. 

1976 

MAIO 

North-South 

Love 

5072 

Aug. 

2. 

1976 

MAIO 

Vertical 

Rayleigh 

5053 

Aug. 

8. 

1976 

MAIO 

Vertical 

Rayleigh 

5100 

Aug. 

8, 

1976 

MAIO 

Vertical 

Rayleigh 

5118 

Aug. 

8, 

1976 

MAIO 

North-South 

Love 

5100 

Aug. 

14, 

1976 

MAIO 

Vertical 

Rayleigh 

5103 

Aug. 

14. 

1976 

MAIO 

North-South 

Love 

5103 

TABLE  14 


Seismogram  Components  Used  for  the 
Path  Tangshan-Taipei 


Earthquake 

Recording 

Station 

Component 

Wave  Type 

July 

28, 

1976 

TATO 

Vertical 

Rayleigh 

July 

28, 

1976 

TATO 

Radial 

Rayleigh 

July 

28, 

1976 

TATO 

North-South 

Love 

July 

30, 

1976 

TATO 

Vertical 

Rayleigh 

July 

30, 

1976 

TATO 

Radial 

Rayleigh 

July 

30, 

1976 

TATO 

North-South 

Love 

Aug . 

1, 

1976 

TATO 

Vertical 

Rayleigh 

Aug . 

1, 

1976 

TATO 

Radial 

Rayleigh 

Aug  . 

1, 

1976 

TATO 

North-South 

Love 

Sept. 

6, 

1976 

TATO 

Vertical 

Rayleigh 

Distance 

(km) 

1661 

1661 

1661 

1670 

1670 

1670 

1656 

1656 

1656 

1672 


TABLE  15 ; Model  Fit  for  Path  Tangshan-Taipei 
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Figure  Captions 


Figure 

1.  Crustal  thickness  and  crustal  plate  subdivision  of  China 

2.  Physical  features  of  China  and  surrounding  regions  (after 
York  et^  aK  , 1976) . 

3.  Surface-wave  paths  (triangles  are  earthquakes,  dots  are  SRO 
stations) . 

4.  Unrotated  seismograms  for  the  path  Tangshan  - Mashad,  August  8, 
1976  (from  top  to  bottom:  vertical,  NS,  and  EW) . 

5.  Rotated  seismograms  for  the  path  Tangshan  - Mashad,  August  8, 
1976. 

6.  Results  of  multiple  filtering  of  Rayleigh  wave  vertical 
component,  Tangshan  - Mashad,  August  8,  1976. 

7.  Results  of  multiple  filtering  of  Rayleigh  wave  radial  com- 
ponent, Tangshan  - Mashad,  August  8,  1976. 

8.  Results  of  multiple  filtering  of  Love  wave,  Tangshan  - Mashad, 
August  8,  1976. 

9.  Group  velocity  of  the  radial  Rayleigh  component  for  five 
events  of  the  path  Szechwan  - Mashad. 

10.  Group  velocity  of  the  vertical  Rayleigh  component  for  five 
events  of  the  path  Szechwan  - Mashad. 

11.  Rayleigh  wave  group  velocity  for  three  events  of  the  path 
Yunan  - Mashad. 

12.  Group  velocity  of  Rayleigh  waves  for  five  events  of  the  path 
Tangshan  - Mashad. 

13.  Love  wave  group  velocity  for  four  events  of  the  path  Tangshan  - 
Mashad. 

14.  Vertical  component  Rayleigh  wave  group  velocity  for  four  events 
of  the  path  Tangshan  - Taipei 

15.  Radial  component  Rayleigh  wave  group  velocity  for  three  events 
of  the  path  Tangshan  - Taipei. 

16.  Love  wave  group  velocity  for  three  events  of  the  path  Tangshan 
Taipei. 


Figure  Captions  (Cont.) 


Figure 

17.  Theoretical  and  observed  group  velocities  for  the  path  Szechwan  - 
Mashad. 

18.  Theoretical  and  observed  group  velocities  for  the  path  Yunan  - 
Mashad. 

19.  Comparison  of  Rayleigh  wave  group  velocities  for  several  path 
across  China  with  average  continental  dispersion  (Ewing  et  al., 
1957). 

20.  Resulting  shear  wave  velocity  model  for  the  path  Szechwan  - Mashad 

21.  Resulting  shear  wave  velocity  model  for  the  path  Yunan  - Mashad. 

22.  Comparison  of  the  Canadian  Shield  model  (Brune  & Dorman,  1963) 

and  the  Gutenberg  model  (Takeuchi  et  al.,  1964)  with  model  SMI 
and  YM1 . ~ 

23.  Theoretical  and  observed  group  velocities  for  the  Tangshan  - 
Mashad  path. 

24.  Resulting  shear  wave  velocity  model  for  the  Tangshan  - Mashad 
path. 

25.  Comparison  of  the  Canadian  Shield  model  (Brune  & Dorman,  1963) 
and  the  Gutenberg  model  (Takeuchi  et  al_. , 1964)  with  model  TM1 

26.  Observed  and  theoretical  group  velocities  for  the  path  Tangshan  - 
Taipei. 

27.  Resulting  shear  wave  velocity  model  for  the  path  Tangshan  - Tepei. 

28.  Comparison  of  Canadian  Shield  model  (Brune  & Dorman,  1963)  and 
the  Gutenberg  model  (Takeuchi  et  aK , 1964)  with  model  TT1 . 
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